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Abstract 

We consider expectation values of local operators in (continuum) integrable models in a 
situation when the mean value is calculated in a single Bethe state with a large number of 
particles. We develop a form factor expansion for the thermodynamic limit of the mean value, 
which applies whenever the distribution of Bethe roots is given by smooth density functions. 
We present three applications of our general result: i) In the framework of integrable Quantum 
Field Theory (IQFT) we present a derivation of the LeClair-Mussardo formula for finite tem- 
perature one-point functions. We also extend the results to boundary operators in Boundary 
Field Theories, ii) We establish the LeClair-Mussardo formula for the non-relativistic ID Bose 
gas in the framework of Algebraic Bethe Ansatz (ABA). This way we obtain an alternative 
derivation of the results of Kormos et. al. for the (temperature dependent) local correlations 
using only the concepts of ABA. iii) In IQFT we consider the long-time limit of one-point 
functions after a certain type of global quench. It is shown that our general results imply 
the integral series found by Fioretti and Mussardo. We also discuss the generalized Eigen- 
state Thermalization Hypothesis in the context of quantum quenches in integrable models. 
It is shown that a single mean value always takes the form of a thermodynamic average in a 
Generalized Gibbs Ensemble, although the relation to the conserved charges is rather indirect. 

1 Introduction and our main result 

One of the central tasks of many-body quantum physics is the calculation of correlation func- 
tions. In addition to the ground state properties it is also important to consider correlations 
calculated in excited states. These quantities are relevant in (at least) two physical situations: 
they describe the finite temperature properties of the system and they can be used to calculate 
observables after a sudden quantum quench [1]. Besides calculating the space and time dependent 
correlations it is also important to consider the one-point functions, ie. the mean values of lo- 
cal operators. Important examples include the (properly regularized) products of field operators, 
which can describe densities of conserved charges or certain amplitudes of inelastic many-body 
processes [2, 3]. 

In this paper we consider one-point functions in one-dimensional integrable models. These 
theories can be solved by the Bethe Ansatz, which means that the spectrum and also the ther- 
modynamic quantities can be determined exactly. Moreover, there are powerful methods available 
to obtain correlation functions. One of the most general methods is the so-called form factor ap- 
proach, which consists of two main steps. The first step is the calculation of the matrix elements 
(form factors) of the local operators in the eigenstate basis of the Hamiltonian, which can then 
be used to obtain the correlation function as the sum of a spectral series. 

In this work we elaborate on the form factor expansion for one-point functions, more specifically 
we consider mean values in highly excited states with a smooth distribution of Bethe roots. We only 
consider theories with purely diagonal scattering, but otherwise the calculations are quite general. 
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We consider both relativistic Quantum Field Theories (QFT's) and non-relativistic many-body 
systems in the framework of Algebraic Bethe Ansatz. There are three main sources of motivation 
for the present work. 

The first motivation comes from integrable QFT, where the form factors are obtained as 
solutions to the so-called form factor bootstrap program [4, 5, 6, 7]. Originally the form factor 
approach was developed to obtain two-point functions at zero temperature, but it was hoped for 
that the form factors could be used to describe finite temperature correlations too [8, 9, 10, 11, 12]. 
The formal Boltzmann averages are plagued with divergences, which arise from singularities of 
the form factors and different contributions to the partition function. It is fairly easy to show that 
the divergences cancel, but the determination of the left-over pieces is in general a demanding 
problem. 

In [9] LeClair and Mussardo proposed compact integral formulas for the finite temperature 
correlations, which involve the (appropriately regularized) zero-temperature form factors and the 
occupation numbers as determined from the Thermodynamic Bethe Ansatz (TBA). In the case of 
two-point functions the results of [9] were questioned in [10, 13]. On the other hand, it was proven 
it [10] that for one-point functions the LeClair-Mussardo (LM) formula is correct for operators 
describing densities of conserved quantities. Moreover, it was shown in [14] using a finite-volume 
regularization that the LM series is valid for arbitrary operators up to the first three non-trivial 
orders. However, the methods of [14] can only be applied order- by-order, and the calculations 
become rather tedious for higher orders, therefore they are not suitable for an all-order-proof of 
the LM formula. 

The second motivation for the present work comes from recently discovered relations between 
integrable QFT and non-relativistic models of Algebraic Bethe Ansatz (ABA). In [15, 16] a non- 
relativistic limit of the sinh-Gordon field theory was applied to obtain physical observables in the 
Lieb-Liniger model [17, 18]. They performed the non-relativistic limit of the LM series for certain 
sinh-Gordon operators to arrive at the (temperature dependent) local correlations in the Bose 
gas. It was pointed out later in [19] that the non-relativistic limit not only applies to the physical 
observables, but also to the individual form factors entering the integral series. This means that 
the final results of [15, 16] can be formulated using only the quantities calculated in the Bose gas, 
without actually referring to the sinh-Gordon model. Evidently, this suggests that the LeClair- 
Mussardo formalism is not at all limited to relativistic field theory and it could be applicable to 
other Bethe Ansatz solvable models too. We stress that although there is a vast literature devoted 
to correlation functions in the Bose gas, the papers [15, 16] were apparently the first ones to 
calculate the many-body local correlations using Bethe Ansatz techniques. 

It was a vital point of [15, 16] to apply the LeClair-Mussardo formalism in the case of an 
arbitrary prescribed particle density, or in other words in the presence of a non-zero chemical 
potential /i. Although this seems to be quite natural, we remark that neither the original derivation 
in [9], nor the partial proof of [14] applies in this case. Indeed, the calculation of [14] heavily relied 
on the fact that if \i = 0, then the particle density p = N/L can be expanded into a low- 
temperature expansion and that p(T = 0) = 0. On the other hand, the success of [15, 16] suggests 
that it should be possible to prove the LM series using a different approach, which should work 
at arbitrary and both in relativistic QFT and in Algebraic Bethe Ansatz. 

The third motivation for the present work is the quench problem of Fioretto and Mussardo [20] , 
where the authors consider the long-time limit of one-point functions after a sudden global quench. 
The resulting integral series is a generalization of the LM formula, where the statistical weight 
functions do not follow from thermodynamic quantities, but from the microscopic amplitudes of 
the quantum quench and a generalized TBA-like dressing procedure. Once again, the compact 
form of the final result suggests that it should be possible to derive it using alternative methods, 
bypassing the combinatorial difficulties encountered in [20]. 
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In this paper we approach the above mentioned problems using a finite volume regularization 
for the mean values. We start from the Boltzmann average for thermal one-point functions in a 
finite volume L, which can be written as 

where the summation runs over a complete set of states. Note that contrary to the infinite volume 
case the above formula is perfectly well-defined in finite volume. For simplicity we assume that 
there is only one particle type in the spectrum; the generalization to arbitrary diagonal scattering 
theories is straightforward. In this case the Bethe equations (the quantization conditions for the 
finite volume states \i) , ) can be written as 



\S{0j-e k ) = l, j = i... N 



Here S(9) is the elastic S-matrix which is a pure phase in this case. The variables 9j are the 
rapidity parameters and one-particle momenta and energies are given by the functions p{6) and 
e{6). The results to be presented below apply both to relativistic and non-relativistic situations. 

In [14] a low-temperature expansion of (1.1) was performed in massive integrable field theory. 
The derivation was built on the fact that if the chemical potential is zero then the contribution of 
the A-particle sector of the Fock space can be estimated as (mLe~ m ^ T ) N , where L is the volume 
and m is the mass of the particle. Therefore, in the regime 

e m/T (1.2) 

the volume can be large enough to replace sums with appropriate integrals, and in the same time 
the relevant contributions in (1.1) would still come from sectors with a low number of particles. 
This was used in [14] to derive a systematic low-temperature expansion and the L — > oo limit was 
taken at the end of the calculations. The same ideas were employed in the papers [21, 22, 23, 24]. 

In the present work we pursue a complementary approach: we consider a very large volume 
such that there are a large number of particles present in the system. In this case it is expected 
that the thermal average will be dominated by those states in which the distribution of Bethe roots 
is described by the infinite volume TBA equations [25, 26, 27]. In fact the TBA yields densities 
which minimize the free energy functional for the partition function, and one can argue that the 
insertion of the local operator does not shift this saddle point solution. Moreover, one can show 
that in the thermodynamic limit it is sufficient to consider only one particular state and evaluate 
the mean value on this state. 

Let us therefore introduce a "thermal state" as 

\n) L = \e 1 ,...,e N ) L L (n\n) L = i, 

where it is understood that the rapidities Oj satisfy the BA equations, the particle number is 
N = [pL] with p being the infinite volume particle density ([x] denotes the integer part of x), and 
the distribution of roots is given by the TBA equations (detailed definitions will be given below) . 
Then the mean value (1.1) is given by the thermodynamic limit 

(o) T = lim (n\o\n) L (i.3) 

N,L— >oo 

In general the matrix elements with a large number of particles are complicated objects, and the 
evaluation of the thermodynamic limit can be a difficult problem [28, 29, 30]. However, in the case 
of mean values the situation is much simpler. In [14] it was shown that a diagonal matrix element 
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is the sum of its properly defined disconnected terms, which can be expressed using form factors 
with a lower number of particles. Then the remaining task is to perform the thermodynamic limit 
for the disconnected terms. 

It is useful to investigate the mean value in a more general situation. Let us consider a Bethe 
state \ip) L in a large volume with N = [pL] particles. We assume that the distribution of Bethe 
roots is given by the functions p(°'(6) and p( h \8), which describe the densities of occupied states 
and holes, respectively. They satisfy the constraint 



{ e) + pth)(e) = — p \o) + J^ —v{e-e>)p^{6') (1.4) 



where <p = — i4klogS(9) is the scattering kernel and p'(9) is the derivative of the one-particle 
momentum. The particle density is then given by the integral 

dd p(°)(0) 

We consider the thermodynamic limit of the mean value: 

(<D)= lim (^m) L (1.5) 

N,L— too 

For the above limit we find our main result: 

w-Ea/I^-*^))^.-.^) (L6) 

Here F® n c are regularized evaluations of the diagonal form factors (to be defined in the main text) 
and the statistical weight functions are given by 

f{6) = p( h )(e) + p(°)(ey (L7) 

The nth term in the integral series arises from the n-particle disconnected terms of the original 
matrix element (1.5). The form factors describe n-particle processes over the Fock vacuum and 
the only effect of the non-trivial background is the appearance of the weight functions f(9). 
Remarkably (1.6) holds both in massive relativistic Field Theories and in the non-relativistic ID 
Bose gas, which is a gapless system. This will be proven in sections 2 and 3, respectively. In the 
latter case (1.6) can be considered as an (improved) I/7 expansion, where 7 is the dimensionless 
effective coupling of the Bose gas. 

In the finite temperature problem one introduces the pseudo-energy function e{9) as 



pW(9) 

where e{9) is determined by the TBA equation 

tiff 



/oo mi 

—^{e-e')\og{i + e-^\ (1.8) 

with e(9) being the one-particle energy. Applying (1.6) in this case one obtains 

^Es/fr-fr (ndJ 0*,...a), (1.9) 
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which is just the LeClair-Mussardo series as originally proposed in [9]. Note that the general 
result (1.6) allows for an arbitrary chemical potential //. As a special case one can take the zero 
temperature limit with a fixed particle density to obtain 

<°w, L fK<*- < lio > 

where B is the Fermi-rapidity. This remarkably simple and intuitive result was written down for 
the first time in [15, 16]; it applies whenever the ground state with a fixed number of particles is 
unique. 

Let us consider the (typically infinite) set of local onserved quantities Qi with i = 1, 2, . . . , oo. 
The first two members of the series can be chosen as the particle number and the energy. If the 
one-particle eigenvalues are given by the functions qi{6) then the macroscopic values of the charges 
are given by the integrals 

Qi = L J d0p(°\6) QiiO) (1.11) 

Assuming that there is a one-to-one correspondence between the root densities and the conserved 
charges (ie. the relation (1.11) can be inverted) the main result (1.6) can be interpreted as a 
"Generalized Eigenstate Thermalization Hypothesis" [31]: the mean value always takes the form 
of a thermodynamic average and the weight functions only depend on the macroscopic value of 
the conserved charges. 

The remainder of the paper will be devoted to the proofs and applications of our general 
formula (1.6). In section 2 we establish (1.6) by performing the thermodynamic limit of the mean 
value in relativistic QFT. In 2.2 we extend the results to boundary operators in Boundary QFT 
and prove the corresponding LeClair-Mussardo formula proposed in [32]. In section 3 we apply 
our formalism to the Lieb-Liniger model in the framework of Algebraic Bethe Ansatz. In 4 we 
explain how to apply our general formula to quench problems. We re-derive the results of Fioretto 
and Mussardo [20] and we also discuss more general quench situations. Section 5 includes our 
conclusions, and a number of technical details are presented in the appendices. 



2 Form factors and expectation values 

Let us consider a massive integrable QFT. The basic object of the theory is the factorized 
S*-matrix [33, 34]. For simplicity we assume that there is only one particle type in the spectrum 
with no internal degrees of freedom. In this case the S'-matrix is a pure phase: 

S(0i - 0j) = e^-W 

One-particle energy and momentum are given by the functions 

e(6) = mcoshO and p(9) = msinh 6 

In infinite volume the Hilbert space is spanned by the asymptotic states 

\0±, .. . ,0n) 

Their energy and momentum (and higher conserved charges) can be calculated additively. 
Consider a local operator 0(x,t). The form factors of O are defined as [5] 

#m|01) • ■ • > On) — ■ ■ ■ , #m • • • i On) (2.1) 
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With the help of the crossing relations 

f m,n(9ii ■ ■ • >0Af|0i> • • -i^n) = f m-i,n+i(9ii ■ ■ ■ >9m-i\6m + «r>0ij ■ ■ • ( 2 - 2 ) 

N fc-1 

-l,JV-l(^li • • • > ^Af-ll^l> • • • ) 0fc~l> 0fc+l • • • ) @n) 

k=l 1=1 
all form factors can be expressed in terms of the elementary form factors 

F^(e 1 ,...,e N ) = {o\o(o,o)\e u ...,e N ) (2.3) 

which satisfy the following equations: 

I. Lorentz transformation: 

F%{0 X + A, 9 2 + A, . . . , 9 N + A) = exp ( So A) F%(9 1 ,9 2 , ...,9 N ) (2.4) 

where so denotes the Lorentz spin of the operator O. 

II. Exchange: 

F N (9i, ■ ■ ■ ,0 k , Ok+i, ■ ■ ■ , On) = 

S(6k — 0fc+i)-Fjv (^l) ■ • ■ > 8k+i,8k, ■ ■ ■ , &n) (2-5) 

III. Cyclic permutation: 

F%(e 1 + 2iTr ) e 2 ,...,e N ) = F%(e 2 ,... ) e N ,e 1 ) (2.6) 

IV. Kinematical singularity 



- i Res F# +2 (0 + m,e',9x,...,9 N )= (l-T\ S(d' - d k )] F J 
e = 6 V fe=i / 



(2.7) 



There is also a further equation related to bound states which we do not need in the sequel. 

In order to derive a regularized form factor expansion for the mean values we consider the 
theory in a finite volume L. Assuming periodic boundary conditions the Bethe equations are 
written as 

e i Pj L JJ s{ q. _ 6k) = 1 ^ j = l... N, 

k^j 

where pj = p(6j). In the logarithmic form: 

Qj = pjL + ~ h) = 2irlj j = l...N, Ij G Z (2.8) 

k+3 

The multi-particle energies and momenta are calculated additively: 

P N = Y j p(Q j ) E N = Y,e{6 ] ) 

3 3 

In the following we denote the finite volume states as 
where it is understood that the rapidities solve eq. (2.8). 
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It is a very important task to determine the finite volume form factors in terms of the infinite 
volume quantities. In the case of relativistic QFT this problem was solved in [35, 14]. In a generic 
case the result reads 

(e[, ...Mo^,..., e N ) L = fSM^-M^-M l 

The quantities and pm are the multi-particle density of states and they are given by 

PJV (0i, ...,0n) = det J ij , J ij = ^ (2.10) 

The error exponent p is universal in the sense that it only depends on the analytic properties of 
the S-matrix and not on the particular form factor in question. Equation (2.9) is very natural in 
the Algebraic Bethe Ansatz approach for models with si (2) or U q {sl{2)) symmetry [36]. In these 
models (2.9) is exact and p^ is just the Gaudin-determinant, which describes the norm of the 
Bethe wave function. 

The relation (2.9) is valid when there are no kinematical poles present due to colliding rapidi- 
ties. It was pointed out in [14] that there are only two situations when rapidities do coincide in 
a finite volume. One possibility is to have zero-momentum particles in parity symmetric states; 
these matrix elements will not be needed here, therefore we omit the details which can be found 
in [14]. The other case is the problem of diagonal matrix elements (expectation values), which 
will be discussed below. Note that the diagonal limit in (2.9) is ill-defined due to the multiple 
kinematical poles prescribed by (2.7). 

In the following we restate the relevant results of [14]. Consider the diagonal form factor of N 
particles in infinite volume: 

+ e 1 ,...,e N + e N \e N , ...,e 1 ) (2.11) 

Here the singularities have been shifted off by the infinitesimal quantities £j. It was proven in 
[14] that there exists a finite limit when all go to zero simultaneously with their ratios fixed. 
Moreover, there are two special evaluation schemes which respect the physical requirement that 
the diagonal form factors should not depend on the order of the rapidities. First of all, one can 
consider the symmetric limit 

Ff Ns (9i,..., N ) = lim F^ N {9 X +e,...,6 N + e\0 N , . . . , 9 X ) (2.12) 

On the other hand, one can also consider the connected part of the diagonal form factor which is 
defined to be the contribution to (2.11) which does not contain any singular factors of the form 
£i/ej and products thereof: 

Ff NiC (6 u ...,9 N ) = (finite part of) F^fa + e 1 ,...,6 N + e N \9 N , ...,9 X ) (2.13) 

The general structure of the singularities in (2.11) was worked out in [14]. For future use we 
recall the relation between the symmetric and connected evaluation schemes. Firs we introduce 
the necessary notations. 

Let us take n vertices labeled by the numbers 1,2, ... ,n and let G be the set of the directed 
graphs Gi with the following properties: 

• Gi is tree-like. 

• For each vertex there is at most one outgoing edge. 
For an edge going from i to j we use the notation Eij. 
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Theorem 2.1 The function F® N (0±, . . . , 8n) can be evaluated as a sum over all graphs in G, 
where the contribution of a graph Gi is given by the following two rules: 

• Let Ai = {a\,a2, ■ ■ ■ ,a m } be the set of vertices from which there are no outgoing edges in 
G{. The form factor associated to G{ is 



For each edge Ej^ the form factor above has to be multiplied by 

vie, - e k ) 



This theorem was proven in [14] using the kinematical singularity axiom (2.7). 

Now we are in a position to state one of the central results of [14], which is the following 
relation between diagonal form factors in infinite and finite volume: 

. . . , e N \o\e x , e N ) L = 1 £ f£ iS ({mW-„({M) (2.15) 

PN{ 1, ■ ■ ■ , N) | 0+ | U{6) _} 

The sum runs over all possible bipartite partitions. Similar expansions are known for mean values 
in the Algebraic Bethe Ansatz literature [37, 38, 39, 18]. In that context the symmetric evaluation 
is called the "irreducible part" of the form factor. More details about the relation of our expansions 
to those in ABA are presented in section 3. 

A very important result of [14] is that (2.15) can be expressed alternatively in terms of the 
connected evaluation of the diagonal form factors. To restate this result we introduce a restricted 
Gaudin determinant as follows. For a given bipartite partition 

{o 1 ,...,e N } = {e + }u{e-} 

\{8 + }\=N-n and |{#_}|=n 
we define the restricted determinant 



p N - n ({0 + }\{e-}) = detJ + (2.16) 

where J+ is the sub-matrix of J corresponding to the particles in the set {0+}. Note that 
p~N- n ({8 + }\{9_}) still contains information about the complementary set of rapidities {#-}, and 
it should not be confused with the density p7v-n({#+}) which depends only on the rapidities {0+}. 
With these notations, the alternative expression for the expectation value reads 

(e,,..., o N \o\o u . . . , e N ) L = 1 < c ({MW-n({MKM) (2.17) 

PN{ 1, ■ ■ ■ , N) |g + } U { 6 ,_} 

The equivalence of (2.15) and (2.17) is proven in Theorem 2 of [14]. We wish to note that the two 
relations are expected to be exact to all orders in 1/L, but there are residual finite size effects of 
order 0{e~^ L ). The general structure of exponential corrections to form factors is not yet known, 
however it is expected that the growth of the coefficients of the 0(e~ fl L ) terms with any p! > p 
is only polynomial in N. Therefore it is safe to neglect them in taking the thermodynamic limit 
of (2.15) and (2.17). 
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2.1 Evaluation of the thermodynamic limit 

We consider the diagonal matrix element 

(0) NL = (e 1 ,...,e N \o\e 1 ,...,e N ) L (2.1s) 

with a large number of particles in a large volume L. The idea is to evaluate expressions (2.15) 
and (2.17) in the L — > oo limit assuming that the distribution of roots is given by the smooth 
functions p(°\9) and p^(9). They describe the density of occupied roots and holes, respectively. 
They satisfy the constraint 

,M(0) + pW(0) = ^l + /~f^-^ (o) (0'), (2-19) 

which follows from the thermodynamic limit of the Bethe equations. The kernel in (2.19) is given 
by ip = —i-fjfi log S(9) and p(9) is the one-particle momentum. The total particle density is given 
by the integral 

N 



P = T = I dBpM(0) 



Let us introduce the quantities 



, ~ Pn(pi, ■■■,9n) 

\{6-}\=n 



Pn{9\, ■ ■■ , 9n) 

-}\=n 



such that 



N N 



(0)n,l = Y. C ^ = Y. D ^ 

n=0 n=0 

It will be shown that both C n ^ and D n ^ possess a well-defined L — > oo limit for arbitrary n. 
This way the expectation value is expressed as 



(<D)=J2Cn = Y,D r , 



n=0 n=0 

where 

C n = lim C n> L D n = lim D UjL 

L— >oo L—>oo 
We assume that performing the summation over n and taking the thermodynamic limit can be 
exchanged. Note that the total number of C n ^i and depends on L, because n < [pL\. However 
this is not a serious problem, because for any n there will be a large enough volume with N > n, 
so that C nj L will appear in the sum; moreover it will have a finite limit as L is sent further to 
infinity. 

The n = term reproduces the infinite volume vacuum expectation value: 

C ,l = D 0)L = (O|0|O) 

In the present approach this vacuum expectation value is interpreted as the maximally discon- 
nected part of the diagonal form factor. 
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In the case of n = 1 the quantities C\ t L and are expressed as a single sum over 9j 6 
{6»i, . . . , 9n}- 

n sr t?o (Q \PN~i{Oi,---,o j ,...,e N \e j ) 



Pn{Qi, ■ ■ ■ , 9n) 

Pn-i (di, ■ ■ ■ ,0j, . . . , On) 
Pn(0i, • • • j On) 



Let us assume that the ratios of the determinants have a well-defined L — > oo limit while keeping 
6j fixed. Then the summation over 9j can be expressed as an integral: 



£ ^ jdd P (°\e)L 



The ratios of the determinants scale as 1/L, therefore C\ l an d ~D\,l behave as O(L ). This is 
expected in order to have a well-defined L — > oo limit. 

The limiting values of the ratios of the determinants can be determined using the techniques 
of [38, 18]. In Appendix A it is shown that 

PN^i{e 1 ,...,e j ,...,e N \e j ) _ 1 



pjv(0i, ...,9 N ) 27rL(p(°)(9 j ) + pWtfj)) 

PN -i{e l ,...,e jl ...,9 N ) i 

PN (9 U ...,9 N ) 2 7 rL(p(o)(9 j ) + p^(9 j )f { ° 3) 



Here 



W (0) = exp {- f ^-f(9')<p(9-9') 
and the weight function f(9) is defined as 

f (S \ = pio] ( d) 

J{> p(°)(0) + p( h )(9) 
Putting everything together one obtains 



d = lim C 1)L = [ p-mF&iO) 
D 1 = lim D 1>L = [ ^-f(9M9)F 2 ° s (9) 



(2.20) 



One can repeat this line of reasoning for the case when a fixed (n > 1) number of rapidities are 
chosen to be in the subset {9-}, which are to be substituted into a diagonal 2n particle form factors. 
It is easy to see that in the thermodynamic limit there will be no complications, in particular the 
application of the Pauli principle for the summation over {9-} only causes corrections of 0(1/L). 
Therefore one can immediately write down the general result 



and 



Dn= h.j d t--- d t (n/^M^jo^---.^) ( 2 - 22 ) 
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The expectation value of O is then given by 



y- 



n 



(2.23) 



/It- -It IImw^) nO^-a) 



n! , 

fe " \i=i 



The two different expressions represent a rearrangement of terms, similar to the equivalence of 
(2.15) and (2.17). We wish to note that (2.23) is the main result of this work. 

Both series in (2.23) are expected to be convergent. We are not able to present a rigorous 
proof, because the behaviour of the form factors is in general not known. However, it is generally 
assumed that their exists a K 6 K + such that 

\F° +m (e' 1 ,...,e' n \e 1 ,...,e m )\<K n+m for e^ei (2.24) 

This implies that both functions F® n s and c are bounded in the sense of (2.24) and the series 
(2.23) are convergent. 

2.2 Boundary operators in Boundary Field Theories 

It is possible to extend the previous calculations to relativistic Boundary QFT's [40, 41]. 
Consider an operator 0{t) localized at the boundary of a half-infinite system. The boundary is 
supposed to be integrable; the scattering off the boundary is described by the elastic reflection 
factor R(9). Asymptotic states can be defined similarly as in the bulk case; the only difference is 
that for incoming (outgoing) states all the rapidities are positive (negative). Form factors of O 
are then defined as 1 

G%{0 U ...,0 N ) = (0\O\9 1 ,...,d N ) 

The analytic properties of the functions (also called boundary form factor equations) were 
derived in [42]. For the sake of brevity we do not cite them here. 

We will be interested in the mean values of the boundary operator O. We start with a finite 
system of size L with two integrable boundaries o, b with elastic reflection factors R a ^{9)\ the 
operator O lives on the boundary a. The Bethe equations read 

e iQf = e ^L Ra ^ j)Rb ^ ej) Yl S (9 k - 9j)S(9 k + 0j) = 1, 9j > 

Mi 



The rapidities are restricted to take only positive values. The iV-particle density of states is 
obtained as 

^(W) = det l 7 B , J* 



dQf 

89 'j 



Analogously to (2.16) we define a restricted density for a given partition {9} = {#+} U {6L} as 
the sub-determinant belonging to the subset {#+}: 

The relation between finite volume and infinite volume form factors was worked out in [43]. In 
the generic case with no coinciding rapidities they take the same form as (2.9) with the obvious 



1 Form factors of boundary operators are denoted by G% to distinguish them from the form factors of bulk 
operators. 
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replacements FS — > and p^ — > pfj- For diagonal form factors one obtains 

(e 1 , . . . , 9 N \o\e 1: . . . , e N ) L = 1 £ G& iC ({MK_ n ({M|{0_», (2-25) 

PN {9 1 ,...,e N ) {9+}u{9 } 

where G^j C are the connected parts of the diagonal form factors defined analogously to (2.13). 
Clearly, (2.25) is a simple generalization of the relation in the bulk (2.17). It is also possible 
to define the symmetric evaluation of the boundary form factors as in (2.12), but the naive 
generalization of (2.15) does not hold in the boundary explained in detail in [43]. 

We are interested in the thermodynamic limit of (2.25) along the lines of the previous sub- 
section. Using the techniques of Appendix A it can shown that the behaviour of the ratios of 
determinants is the same as in the bulk. For example 

lim p^(fli,-.-,flj,...,0Ar|flj) 1 



where p(°) (9) and p^ (9) describe the densities of occupied roots and holes, respectively. Therefore, 
the thermodynamic limit of (2.25) results in the series 



(O) 

where 



§-f £ (n /<*>)<*,<*.-.«.>. e-»> 



p(°)(9) + p( h ){9) 

In the case of finite temperature expectation values one can show that the distribution of 
roots will be given by the same functions as in the periodic expected by thermodynamic 

arguments. Therefore 

(°h, = Ehr^--- I ^inr^rK^-.U (2-27) 

where e{9) is the solution of the TBA with periodic boundary conditions (1.8). The series (2.27) 
was proposed in [32] and it was proven there up to third non-trivial order using the techniques of 
[14]. The all-orders proof together with the general formula (2.26) is a new result of this work. 




3 Application to the ID Bose gas 

In this section we consider the ID interacting Bose gas also known as the Lieb-Liniger (LL) 
model [17, 18]. The second quantized form of the Hamiltonian in volume L with periodic boundary 
conditions is given by 

#ll = J dx (d x tfd x V + 0*%%*) . (3.1) 

Here *f?(x, t) and ^Sf'{x,t) are canonical non-relativistic Bose fields satisfying 

[*(x,t),9\y,t)]=6(x-y), (3.2) 

and c is the coupling constant. We used the conventions m = 1/2 and H = 1. The Fock vacuum is 
defined as 

*|0)=0, (01^=0 
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The eigenstates of the Hamiltonian (3.1) can be constructed using the Bethe Ansatz. The scat- 
tering of the particles is described by the two-particle S-matrix 

X-ic 
A + ic 

where A is the non-relativistic rapidity variable and multi-particle energies and momenta are given 
by 

e n = ^2 A i Pn = a j 

3 3 

The Algebraic Bethe Ansatz (ABA) provides a framework to calculate form factors and corre- 
lation functions. However, before turning to the ABA solution we present the results of the papers 
[15, 16], where a non-relativistic limit of the sinh-Gordon theory was performed to obtain physical 
quantities in the LL model. They considered the temperature dependent expectation values 

(0 fc > = /*tV\=0 fc (T,c) p k 



where p = N/L is the particle density. The dimensionless quantities are important for the 
phenomenology of the Bose gas, for example g% describes the recombination rate of the gas [2, 3]. 
The main result of [15, 16] is the integral series 



E l f dXi dX N ( -j^r 1_ 

N=k \J=1 



N= 

where e(A) is the solution of the non-relativistic TBA equation 

^(e-e^logil + e-^), (3.4) 

where 

^( A ) = log 5*LL 



8X & ULJ c 2 + A 2 

The functions F^ N c appearing in (3.3) were derived by the non-relativistic limit of certain form 
factors in the sinh-Gordon model. The first few examples were given explicitly in [16]. They vanish 
for N < k and the asymptotic behaviour at c — > oo is given by 

In the following we re-establish the series (3.3) using the methods of Algebraic Bethe Ansatz 
[37, 38, 39, 18]. Also, we clarify the connection between the functions F^ N c above and the form 
factors calculated from ABA. 

In the ABA approach the central object is the monodromy matrix, which can be written as 

The commutation relations satisfied by the entries can be expressed in a compact form as 

R(X, n)T(X) <g> T{n) = T(n) <g> T(X)R{X, //) (3.6) 
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Here R(X, //) the R-matrix is of XXX-type: 



(f{\,n) \ 

g(\,n) 1 

1 ff(A,/i) 

V /(A,/i)/ 



where 

/(A, At) = — : ff(A, M) 



X — fi A — /i 

The vacuum eigenvalues of the operators -A(A), C(A) and D(X) are 

C(A)|0)=0 A(A)|0) =a(A)|0) D(X)\0) = d(A)|0) 

For definiteness we specify 

a(A) = e~ iAi/2 d(A) = e iAL/2 

The function /(A) = a(A)/d(A) = e~ tXL will be used extensively. 
The normalized operators B(A) and C(A) are given by 

B(A) = B(X)/d(X) C(X) = C(X)/d(X) 

The Bethe states are then defined as 

N 

|A 1 ,...,A A r} L = n B (^)l°) ( 3 - 7 ) 

These states are eigenvectors of the Hamiltonian (3.1) if the rapidities satisfy the Bethe equations: 

l(Xj) J] S LL (X k - A,) = 1 (3.8) 

Mi 

For later use we also introduce the logarithm of the l.h.s. as 

Qj(X) = ilogZ(A) +iJ^logS LL (A - A fe ) = AL + 2^atan A °_ - 

kj^j k^j 3 k 

We will be interested in the matrix elements of the operators 

O k = (^) k ^i k 

First of all let us define the un-normalized "off-shell" matrix element as 

N N 

M% N {{X C } N , {X B } N , {l c } N , {l B } N ) = <0| [J C(Af )(*t)*** J] B(Af )|0) (3.9) 

j=l k=l 

Here it is understood that the r.h.s. is calculated using the commutation relations (3.6) only; the 
Bethe equations (3.8) are not assumed. This way the function Mjy depends on AN variables given 
by Xf, Xf and if = l(Xf), if = l(Xf). Note that M k N = 0ifk>N. 

The "on-shell" matrix element or "form factor" is defined as a special case of (3.9) when the 
two sets of rapidities are solutions to the Bethe equations. In this case the variables lj can be 
expressed in terms of the functions /(Aj, Xj) and one obtains the function 

^2n({^ C }n\{X b } n ) = M k N ({X c } N , {X b }n, {1 C }n, {1 B }n) 
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which only depends on 2N variables. 

It can be shown that the form factor has first order poles in the case of colliding rapidities. 
The residue at A^ — > is given by [39, 44] 



^2N^ C }n\{X B }n)\ X C^ x 



IC 



N-l N-l 



N A N 



I fjxfXi ~ I f'xfx, I >< ^2N-2({X C }N-l\{X B }N-l) 



3=1 3=1 



Here f.^ = f[Xj ' — X N ' ). Other cases A^ — > \j follow from the symmetry properties of the 
form factor. Note that equation (3.10) does not depend on k, ie. the singularity properties are the 
same for all Ok- Moreover, it can be shown that an analogous equation holds for the operators 

Ojfe,i = (* t )** 1 , MeN 

In Appendix B we present a general proof of (3.10). 

Note that (3.10) has essentially the same structure as the kinematic singularity axiom (2.7) 
in relativistic QFT. There are two differences which are related to different normalizations of the 
Bethe vectors; this was pointed out in [19]. First of all, the norm of the states (3.7) is equal to 



(Ai, . . . , Ajv|Ai, . . . , X N ) = c N [T[fj k )p N (\ 1 ,...,\ N ) (3.11) 
Here /0jv(Ai, . . . , An) is the Gaudin determinant: 



pjv(Ai, • . • , Aat) = det J, J jk 



On the other hand, in relativistic situation the norm of the finite volume states is given simply 
by the density of states pjy, as it is implicitly assumed in (2.9). 

The other issue is related to the exchange of rapidities in the Bethe state. The B(X) operators 
commute with each other, therefore the states (3.7) are totally symmetric with respect to the 
exchange of rapidities. This property also applies to the finite volume matrix elements in relativistic 
QFT, as it was emphasized in [35]. On the other hand, in relativistic QFT the infinite volume 
form factors satisfy the exchange property (2.5), which follows from the Faddeev-Zamolodchikov 
algebra. In the present case the S-matrix is a pure phase, therefore one can introduce the phase 
structure by multiplying with factors of yj 5(A). With a slight abuse of notation the differences 
between the normalization conventions can be summarized as [19] 

N 

|Ai, . . . , Xn)aba ~ ( II ^' fc ) l^ 1 ' '"' ^ N ^QFT (3.12) 

j<k 

In order to make contact with the form factors appearing in the series (3.3) we define 

N _ x 

F 2 k N ({X C }N\{X B }N) = (c N J] ffkffk) ^2n{{X C }n\{\ B }n) (3.13) 

j<k 

The functions F^ N satisfy the exchange property (2.5) and the non-relativistic version of the 
kinematical singularity axiom (2.7). 

Now we are in a position to define the diagonal limits of the non-relativistic form factors. In 
accordance with (2.12) and (2.13) we introduce 

F^Ai, . . . , Ajv) = lim ^(Ai +£,..., Ajv + e\X N , Ai) (3.14) 
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F 2N,c(^h • • • , Ajv) = (finite part of) F$ N (\i +e 1 ,...,X N + e N \X N , . . . , Ai) (3.15) 

The properties of these functions are the same as in the relativistic case. Both are completely 
symmetric in their variables and the relation between F£ N and F% N is given by Theorem 2.1. 
We remark that the normalization in (3.13) commutes with the diagonal limit, because the extra 
factors do not introduce any new poles. 

Now we turn our attention to the mean values of the operators Ok ■ One can define the diagonal 
limit of (3.9) as 

M^(Wn, {1}n,U}n) = lim M&({\ c }n, {X B }n, {1 C }n, {1 B }n) (3.16) 
Formally it will depend on 3iV independent variables {A}at, {1}n and {z}n, where 

z(A) = i^logZ(A) 

The expectation value is then given by the special case of (3.16) when the rapidities satisfy the 
Bet he equations; the remaining task is to express it in terms of the form factors defined in (3.14) 
and (3.15). The result for the normalized expectation value is given by the following formula: 

(O k ) N = , 1 > £ F 2 ^ s ({X^) PN . n ({\ + }) 

(3.17) 



£ Pj lfC ({A_})^_ B ({A + }|{A_}) 



PN (x l ,...,x N ) {A+}U{A _ } 

Clearly, the above formulas are the non-relativistic versions of (2.15) and (2.17). The equivalence 
between the first and the second line can be proven with the non-relativistic version of Theorem 
2 in [14]. Although (3.17) could be proven by performing a non-relativistic limit along the lines 
of [15, 16], we felt it worthwhile to provide a derivation using only the techniques of ABA. This 
is presented in Appendix C. 

Performing the thermodynamic limit of (3.17) along the lines of subsection (2.1) one readily 
derives the integral series (3.3) and also the alternative form 

« - E m / 1 ■ ■ ■ t (n rr^r) F - (A ■ x »> (3 ' 18) 

where 



N " \j=l 



.(A)=e X p(-/f TT l wV (A-V)) 



Equation (3.18) is a new result of this work. It is reminiscent of the integral series for the two-point 
functions derived in [37, 38, 39]. In fact these works consider the thermodynamic expectation value 

(e aQ ^ , a G C, 

where 



Q(x) = / &(x)V(x) 
Jo 

Expectation values of (powers of) Q(x) are then obtained as 

= Q 2 (x). (3.19) 



_^_ e aQ(x) 
da 



ff 2 



Q(x), T^e a «W 



a=0 
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The two-point functions are obtained further by differentiating w.r.t. x. For example the current- 
current correlator is given by 

(J(x)J(O)) = -i^Q 2 (x), J(x) = tf(xMx) (3.20) 

In [37, 38] a form factor expansion was derived for the above correlator. The resulting series has 
similar properties to our formula (3.18), in particular the weight functions are exactly the same. 
On the other hand, the two-point function is a more complicated object because its x-dependence, 
and the integral series in [37, 38] involve a non-trivial dressing procedure for the momenta of the 
multi-particle excitations. 

It seems natural that our series (3.18) could be obtained by taking the (properly regularized) 
x — > limit of the expectation values of Q k (x). In Appendix E we show how this limit works for 
the finite volume matrix elements (3.17) in the case of k = 2 and N = 2. In principle the x — > 
limit could be performed for arbitrary k and N, then the thermodynamic limit would follow in a 
straightforward way. However, we do not pursue this problem here, because the intention of this 
work was to derive the results for the one-point functions in a direct way. 

In order to evaluate the series (3.3)-(3.18) one has to determine the form factors F^ N c and 
F% N ■ The first few cases were derived in [15, 16] using the non-relativistic limit of the sinh-Gordon 
form factors. In principle this procedure can be performed for any k and N, moreover it is easy to 
implement it with symbolic manipulation programs. It was demonstrated in [15, 16] that in the 
strong coupling regime it is sufficient to consider only the first few terms in the integral series. 
However the general forms of F£ N and F% N are not know, therefore at present the exact result 
(3.3)-(3.18) can be considered a formal expansion. 

In Appendix D we derive the form factors for k = 1, 2 and arbitrary N; this gives an indepen- 
tent confirmation of the results of [15, 16] in the cases k = 1, 2 and N = 1, 2, 3. On the other hand, 
the general results are not known for k > 3. One possibility is to solve the recursive equations 
(3.10) or to employ the techniques of [45, 46] to obtain determinant formulas for the off-diagonal 
form factors; the diagonal limit should be taken afterwards. However, this problem is beyond the 
scope of the present work. 



4 Application to quench problems 

In [20] the authors considered the real-time evolution of the expectation value of a local 
operator in integrable QFT after a certain type of quench, which changes the Hamiltonian from 
Hq to H, where H (possibly also Hq) is considered to be integrable. The main assumption is 
that the initial state of the system (which is the ground state of Hq) can be expanded in the 
multi-particle basis of the integrable Hamiltonian H in the form of a boundary state [41]. In the 
simplest case with only one particle type in the spectrum the corresponding expression is 



|B) = exp (| ^-K(d)A(-0)M0))\0) 



where K(6) is an arbitrary function satisfying K{9) = S{26)K{—9). The time evolution of an 
expectation value is then given by 

(0(0, t)) = {B\e iHt O(0,0)e~ iHt \B) 

The main idea is that in the t — > oo limit (and taking a proper time average) only the diagonal 
matrix elements contribute and the off-diagonal ones can be neglected. The authors then arrive 
at the final expression 
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where the weight function is given by the solution of the TBA-like equation 

- log \G(6)\ 2 = - log \G (6)\ 2 - J — ^(0 - 9') log (1 + \G(8')\ 2 ) (4.2) 
The source term is given by 

G = e- 2mTocoshe K(9), 

where tq is a UV cut-off, which can be interpreted as the extrapolation length known in the theory 
of critical boundary phenomena [47]. From a formal point of view a finite To is needed to have a 
normalizable boundary state. 

Here we re-derive equation (4.1) from our general result (1.6). It will be shown that (1.6) is 
applicable with the weight function 

m _ i<wi' 



l + IGW 

First we consider a large volume L. In finite volume the boundary state is given by [48] 

oo 

\b) l = J2 2 M'(Oi,---,0N)K(e 1 )...K(e N )\-e 1 ,e 1 ,...,-e N ,e N ) L 

N=oe 1: ...,e N 

Here the summation runs over all parity symmetric configurations, where the set {— 8\, 9\, . . . , — 9n, 9n} 
is assumed to satisfy the 2iV-particle Bethe equations and we require 0j > 9j for i > j. The ad- 
ditional normalization factors J\f(9±, . . . , 9n) are specific to the finite volume situation. They are 
given by 

Kr(a a \ V P2n(-0i,0i, . . . , -9n, On) 

M . . . , 9 N ) = r— j- T-r , 

PN\yi, . . . , Vn) 

where p2N is the usual 2A r -particle density of states, and p~N is the constrained density in the 
space of rapidity pairs [48]. It can be shown that in the thermodynamic limit J\f will converge to 
a finite value [49], and it will drop out from the calculation of the expectation value. 

The boundary state is normalizable with any finite To and the (time averaged) expectation 
value is given by 

(O) = i, m an I f T « ^f^mf^k ,, 3) 

We assume that in the large time limit it is sufficient to consider only the diagonal matrix elements. 
In this case the calculation boils down to a simple statistical average: 

(0> h aKfff t (4 . 4) 

where for simplicity we denoted 

\i) L = |-0i, 0i. • • • > -On, 0n) l wi = AT(0i, . . . , 9 N )K(9 1 ) . . . K{9 N ) 

This statistical average (4.4) can be considered as a grand-canonical ensemble for particle pairs 
(—0,0). In fact, one can identify the weights as 

|G o (0)| 2 ~e- £(e) 

where we set T = 1 and E{9) is interpreted as a bare excitation energy: 

E{9) = AmT cosh 9- log \K{9)\ 2 (4.5) 
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For any regular K(6) this bare energy has the typical features of a kinetic term: it is a regular 
function for a small 9, whereas for large \6\ it goes to infinity. 

One can now apply standard arguments to look for those configurations which dominate the 
average (4.4). Repeating all the steps which lead to the Boundary TBA equations [50] it is then 
a standard exercise to arrive at (4.2). In fact, (4.2) is identical to the BTBA equation of [50] with 
the substitution R = tq. 

With this we have shown that the distribution of roots is governed by (4.2), therefore (4.1) 
follows from our general formula (1.6). Let us make some further remarks about this result. 

The most remarkable property of (4.1) was already pointed out in [20]: although the boundary 
states only include rapidity pairs, the form factors in the final result represent processes with single 
particle excitations. This can be understood by recalling how we derived our main formulas in 
subsection 2.1. There it was shown that the individual contributions of the LeClair-Mussardo 
formula arise from disconnected terms of a matrix element with a large number of particles. 
Moreover, the thermodynamic limit of the mean value does not depend on the details of the finite 
volume state, but only on the distribution of roots. Therefore it is completely irrelevant that the 
actual state only involves particle-pairs. The distribution of roots is always parity symmetric, 
therefore the final result is completely consistent with the structure of the boundary state. 

To conclude this section we discuss the implications of the main result (2.23) for more gen- 
eral quench situations. Consider the time-evolution generated by an integrable Hamiltonian H 
from an initial state \&o, which does not have the form of a boundary state. Neglecting possible 
complications due to degeneracies, the infinite-time limit of local observables will be given by the 
Diagonal Ensemble 

(O) = ^ffi , (4.6) 
where the summation runs over the eigenstates of H and 

c a = (a|* ) 

Typically the overlaps are not known, therefore it is very hard to determine which states will have 
an important effect on the average (4.6). On the other hand, in the thermodynamic limit it is 
natural to assume that the dominant states will have a smooth distribution of roots. 

Consider such a state \a) and the mean value {a\0\a) L . According to (2.23) it always takes the 
form of a thermal average with the weight functions determined by the root densities. Assuming 
that there is a one-to-one correspondence between the distributions and the infinite set of conserved 
charges this means that the mean value only depends on the macroscopic value of the conserved 
quantities and not on the details of the state a. This result is a generalized form of the "Eigenstate 
Thermalization Hypothesis" [51, 52] appropriate to integrable theories [31]. Assuming further that 
the dominant states in (4.6) will have the same set of macroscopic conserved charges we conclude 
that the average (4.6) can be substituted by a single mean value {a\0\a) L and then the integral 
series (2.23) applies. This result can be interpreted as a Generalized Gibbs Ensemble (GGE) as 
proposed in [53], although the relation to the conserved charges is rather indirect. 

In order to apply these results in other quench situations the assumptions about the root 
distributions have to be justified. More specifically it needs to be checked whether all relevant 
states in (4.6) can be described by a single distribution. This can be a quite challenging problem. 
We wish to stress that our approach only applies to the thermodynamic limit with the prescription 
(4.3) and it is not clear how to obtain finite size effects or observables at large but not infinite 
times with the present methods 2 . 

2 In the problem of Fioretto and Mussardo expectation values at finite times can be obtained using the methods 
of [48]. However, so far only the first few terms have been calculated explicitly and the general all-orders result is 
not yet available. 
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5 Conclusions 



In this paper we studied mean values of local operators in Bethe states with a large number of 
particles and a smooth distribution of roots. The main result is formula (2.23); it can be considered 
as a generalization of the LeClair-Mussardo formula originally proposed in [9]. We showed that it 
applies both to relativistic field theory and the non-relativistic Lieb-Liniger model. The individual 
terms in the series arise from disconnected pieces of the diagonal matrix element; the mth term 
represents m particle processes over the Fock-vacuum. 

Our results are analogous to the expansion for the two-point functions in the Bose gas derived 
in the papers [37, 38, 39]. The finite volume formulas (3.17) are a new result of this work and 
they serve as a starting point for the thermodynamic limit. The functions s in the first line of 
(3.17) (the symmetric evaluations of the diagonal form factors) are analogous to the "irreducible 
parts" of the correlators known in the ABA literature. On the other hand, the expansion in the 
second line of (3.17) employs the functions F' m( . which are the connected parts of the infinite 
volume diagonal form factors. Relativistic Field Theory motivated the use of the connected parts 
and they seem to be new in the context of Algebraic Bethe Ansatz. 

It would be interesting to extend the present approach to non-diagonal scattering theories, 
which correspond to nested Bethe Ansatz systems. In integrable Field Theory the infinite volume 
form factors are explicitly known in a number of models [5, 54], but the relation to the finite volume 
form factors has not yet been determined. The first step would be to write down the corresponding 
generalization of (2.15) for the finite volume mean values. If such an expression is found, then 
the thermodynamic limit could yield a LeClair-Mussardo formula for non-diagonal scattering 
theories. On the other hand, it is not clear how much can be achieved in non-relativistic (nested) 
Bethe Ansatz systems. At present the only available results concern the spectrum (including the 
thermodynamics [55]) or the norms of Bethe wave functions [56, 57], and it is not clear how to 
obtain the form factors. 

In Section 4 we showed that our general result also applies to quantum quench problems and 
we provided an independent and derivation of the results of Fioretto and Mussardo [20]. In the 
context of quench problems our expansion (2.23) can be interpreted as a proof of a "Generalized 
Eigenstate Thermalization Hypothesis". It would be interesting to check whether the assumption 
of the single dominant root distribution holds in other quench situations. 

In this work we only considered the one-point functions of local operators. It is very natural 
to suspect that similar results hold for multi-point correlation functions, however our methods do 
not apply in that case. Concerning the non-relativistic Bose gas it can be shown along the lines of 
[37, 38, 39] that equal-time two-point functions in excited states only depend on the distribution 
of roots, however there is a non-trivial dressing involved for the intermediate state momenta. It is 
not clear how to obtain similar results in the relativistic setting. In the case of finite temperature 
two-point functions (with zero chemical potential) the first few terms of the form factor expansion 
have been derived recently in [24], but it is not yet clear whether the series can be re-summed to 
a compact expression. 
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Appendix 



A Ratios of Gaudin determinants in the thermodynamic limit 

In this appendix we consider the thermodynamic limit of the ratios 

Pn-i(0i, . . . ,9j, . . . , 9n\9j) pn-i{6i, . . . , flj, . . . , On) 

Pn(0±, ■ ■ ■ ,0n) Pn(@i, ■ ■ ■ ,0n) 

The determinants above are defined in the main text. 

As a first step we write the matrix J defined in (2.10) as 



J = G@ where 



_ <p(8jk) 

l 3 ft . 



'3 

N 



■dj = mL cosh 0j + ¥>(#;' fc) 



i=l 



With this notation 
and also 



Pn{Qi-, • • • i On) = det Gn det G 



N 



pN-i(0\, . . . ,9j, . . . ,0 N ) = detGjv-i dete^-i 
The matrix is diagonal, therefore the factorization property also applies to the sub-determinant: 

pN-i(0i, ■■■ ,0j, ■■■ ,0 N \0j) det G N -i det 6at_i 

where Gn-i is obtained from Gn—i simply by erasing the row and column corresponding to 0j. 
First we consider the ratio 

Pn-i (01) • • • , 0j, . . . , 0n\0j) det ®n-i det Gn-i 1 detGAr_i 
Pn{0\, . . . ,0n) detGjv det Gat det Gat 

In the L — > oo limit we have from (2.19) 

•&j -> 2irLp(Q j ) 
The elements o{ Gn can be written asymptotically as 

_ 1 f(0jk) 

ij " /• mop 

The limit of det Gat is given by the Fredholm determinant 

det fi - —K 
V 2vr 

where 

(K(f))(x) = J ^-y( X -y)f{y) 

Intuitively it is clear that this determinant should not change if we erase one "discretization point" 
given by Qj. This would then result in 
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However, we can be more precise about this. Let us decompose detGjy as a single sum of sub- 
determinants along the column j: 



( 1 \ N 1 

detG N = (l - rtO)^) det^-! - E^iO^y detG 



A 7 " 

(0 

iV-1 



There are O(L) terms in the sum on the right hand side, however the sum itself is still of 0(1/ L) 

because the sub-determinants detGjy_ 1 are 0(1/ L) too. This proves (A.l) and for the ratio in 
question we obtain 

pN-i(0i, . . . , 0j, . . . , 0jsi\6j) _ 1 

im PN (e 1 ,...,e N ) ~ 2irL P (e j ) 

The ratio 

pN-i{®i-> ■ ■ ■ ,0j,... ,8n) _ det0Ar_i detGAr_i 

Pn(&i, ■ ■ ■ i &n) detGjv detG/v 

is more involved because 0jv-i and Gn-i are not restrictions of Gat and Gn; rather we have 
to subtract all contributions associated to 9j. Following the corresponding calculation in [38] we 
write 

N #l 

dete y = tf i dete jy - 1 n tfi _ y(git) 



In the thermodynamic limit we have 



det e*_ a _ 1 " / tpiOfl) \ i , , a , 



detO N Lp(6 j )f = \\ Lp(9 l )J 27rLp(6 j ) 



where we introduced the function 



oj(0) = exp ^— J 



rift' 



It is not hard to convince ourselves that 



det Gat \L 
Putting everything together 



detG "- 1 1 + (A.2) 



lim 



p N (9 1 ,...,e N ) 2irLp(6j) 

B Form factors in the ID Bose gas 

Here we present a detailed derivation of the singularity property (3.10). Consider the scalar 
product of two arbitrary states (not necessarily eigenvectors): 

N N 

(0\l[C(\f)l[M(\Z)\0) 

j=l k=l 
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The scalar product has a simple pole as A^ — > A^: 

N N 

k=l 



3=1 
ic 



N-l 



N-l 



N-l 



(B.l) 



N ~ A N 



(In -In)\.1 fNjf^j \x(0\U C(Af ) J] B(Af )|0> 



mod 



i=i 



fc=i 



Here the scalar product on the r.h.s. has to be calculated with the modified vacuum eigenvalues 

a mo d(A) = a(A)/(A, X N ) d mo d(A) = d(X)f(X N , A) 

It should be noted that the remaining rapidities {Ai, . . . , A^v-i} satisfy the modified Bethe equa- 
tions with l mod = a mod /d mod . 

Equation (B.l) is valid for arbitrary vacuum eigenvalues a(A) and d(A). In the physical case 
the residue vanishes because 



Tc 7b( 1 N ~ 1 n) -> c l NZN 

A N ~ N 



where we defined 



z(X) bg Z(A) 



In the Bose gas one has z(X) = L, however it is useful to leave z(X) unspecified. 

In the diagonal case, when X^ — > X^ for every j the scalar product (the norm) depends 
explicitly on the variables lj = l(Xj) and Zj = z(Xj). The dependence on zn is linear and one 
obtains from (B.l) 



d 

dz N 



N N 

<oi n c ( A f ) ii b ( a 



N-l 



N-l 



N-l 



cIn\U x (0| n C(Af ) [J B(Af )|0) mod (B.2) 



k=l 



.3=1 



3=1 



k=l 



Equation (B.2) was used in [58] to prove the norm formula (3.11). 

Let us consider the action of the field operator ^ = ^(0) on Bethe states. It is given by [39] 



N 



N 



( 



k=l 



*UB(Afc)|0) = -»y^Z(A* 
k=i 

We also need the action of 



JV 



[ f(Xk,Xr 



. m=l 

\m^k 



N 



I B(A r 



(B.3) 



N 



( 



(0|IJC(A fc )*t=iv^X; 



k=l 



k=l 



N 



\ 



I f(X m ,Xk) 



. m=l 

\mj^k 



m=l 
J mj^k 



N 



<0| ] [ C(A r 



(B.4) 



m=l 
/ mj^k 



Note that contrary to (B.3) the function l(X) is not present in the pre-factor in (B.4). This is due 
to the normalization of the operators B(A) and C(A). 

We also need the multiple action of the field operators, which is given by 



N 



m N—m 



N-m 



4*II B (Ai)|0> = -H)c U l ( X P ] [ f( X l'K) ; [ K(A-)|0> (B.5) 



3=1 



{A+}U{A-} \j=l 



K j=l 0=1 



0=1 
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Here the summation runs over all bipartite partitions {\}n = {A + } m U {A }N-m and we made 
use of the identity 



ml 



a \i<j 



where the summation runs over all permutations a G S m . The generalization of (B.4) follows 
straightforwardly. 

The analytic properties of the (un-normalized) matrix elements 

N N 

M k N ({X c } N , {X b }n, {1°}n, {1 B }n) = (0| J] C(Af )(*t)*** J] B (Af )|0> 

j=i k=i 

follow from (B.5) and the basic formula (B.l). The residue at A^ — > X B is given by 
M k N ({X°}N, {\ B }n, {1°}n, {1 B }n)\ x c^ x b -+ 

\C _ \B ( l N ~ 1 n) \ II fNjfNj] X M N~l({X C }N-l, {\ B }n~1, {'SodliV-l, {'modl^-l) 
A N A N \j=l I 

The modified /-functions are given by 

^mod(A) = iWjf^ (B.7) 

The form factors F^r are obtained from M% in the case when two vectors are solutions to the Bethe 
equations and the variables lj are expressed in terms of the functions /(Aj,Aj). The singularity 
properties then follows from (B.6): 



N-l JV-1 



¥ n({ a }N,{X }n)\\C^ x b 



IC 



\\ fjNfNj ~ 11 f?NfNj I x ¥ N-l({ xC }N-l, {A B }at_i) 



A N A N \ j= i j= i 



Here we relied on the fact that the rules for the modification of I - ' on the r.h.s. of (B.6) auto- 
matically produce the N — 1 particle form factor. 

C Mean values in the ID Bose gas 

Here we consider the diagonal limit 

M n ({\} n ,{1} n ,{z}n)= lim M* N ({\ c } N ,{\ B } N ,{l c } N ,{l B } N ), 

which depends on 3N independent variables {A}tv, {1}n and {z}]y. The dependence on z^ is 
linear and it is given by the residue (B.6): 

9 , ffc 



M n ({\} n ,{1} n ,{z}n) 



dz N 

dN I [ fNjfNj I x -^^-i({A}aT-1, {Z m od}AT-l, {^mod}jV-l) 
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The modification rule for the variables Zj follows from (B.7) and is given by 

2mod(A) = z{\) + <p(\ - X N ) (C.2) 

The "on-shell" mean value 

(O k ) N ({\} N , {z} N ) = M N ({X} N , {1} N , {z} N ) 

is a function of 2N parameters {A}jv arid {z}m- The dependence on zm is linear and the coefficient 
is 

d ( N ~ X \ 

(O k ) N ({\} N , {z} N ) = c [ f jN f Nj x (O fc )iv-i ({A}at-i, {z mod } N -i) (C3) 



dz N 



The rule for the modification of the variables Zj is given by (C.2). 

It is convenient to define the symmetric evaluation of the diagonal form factor as 



F&v .(W") = lim F*v({A B + £ } N , {X B } 



N, 



The function s is the "irreducible part" of the diagonal matrix element; it differs from the 
defined in the main text only in the over-all normalization. 

Theorem C.l The irreducible part is equal to the mean value at Zj = 0: 

V k 2 N , s ({\}N) = (O k ) N ({\} N , {0} N ) (C.4) 

Proof It is clear from (B.6) that the zjy dependence of the mean value arises from the rapidity 
dependence of the function /(A). Therefore, the irreducible part can be obtained by putting lj = ij 
where {l B } solves the Bethe equations: 

(O k ) N ({\} N , {0} N ) = \imM$ N ({\ B + e} N , {X B } N , {l B } N , {l B } N ) 

£->0 

On the other hand, the form factor is obtained by 

*W{A}at) = lirr iMl N ({X B + e} N , {X B } N , {l B } N , {l B } N ), (C.5) 

where the Bethe equations are satisfied. This means that {I b }n and {I b }n will be given by the 
appropriate products of S'-matrices. However, a constant shift in the rapidities yields 

[P}n = {1 B }n 

This proves the theorem. ■ 
For any partition {A} at = {A + } n U {\~}n-u we define the function 
SV({A + }n, {\~}n-u, {z}n) = 

(n N-n \ / \ 

n n i/( a , + - x ~k )i 2 n i w - a * )i 2 {*> 
j=l fc=l / \l<j<k<n J 

The dependence on zjy is given by 
d 



dz N 



SW({A + } n , {A }n-u,{ z }n) 



N—i f , f \ v q„ t (f\-\„ /, ,\ ,A ;f i„cn-l„ (C- 6 ) 



\ C \Ilj=l fjNfNj) X S N - n ({\ }N-n,{ z mod}N) if A7V G {A } 
"l0 ii\ N £{\+} n 
Now we are in a position to prove the main theorem of this appendix: 



N-n 
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Theorem C.2 The mean value can be expressed as 



(Ok) N = £ F£ iS ({A + }) S N - n ({\-}) (C.7) 

{A+}U{A-} 

where the summation is over the bipartite partitions of the rapidities into two subsets {\}n = 

{\+} n u{\-} N - n . 

Proof The proof is given by induction over N. The form factor vanishes for N < k, therefore the 
first member will be at N = h. In this case there is no dependence on any of the Zj, therefore 
(C.7) is satisfied. Now let us assume that (C.7) is proven for N < M. We prove that it is valid for 
N = M. 

We investigate the Zj dependence of (C.7). According to the assumption of induction it follows 
from (C.3) and (C.6) that the Zj dependence of the l.h.s. and r.h.s. in (C.7) coincides. Therefore 
we consider the value at Zj = 0. At this point the norm functions the r.h.s. of (C.7) vanish and 
one is left with 

KAW) ( c - 8 ) 

It was proven in Theorem 1 that this form factor coincides with the irreducible part of the mean 
value. Thus (C.7) is proven. ■ 
The normalized expression of (3.17) follows from (C.7) by dividing with the norm (3.11). 

D The form factors F^ Ns and F^ N s 

In this section we derive explicit expressions for the functions F^ N and FJ? N for N £ N and 
k = 1,2. 

In the first case the operator in question is 0\ = \E f t\I' which describes the total particle density 
Therefore its mean value is given by 

*t*\ =K= 1 ^» A » (Dl) 
In L Pn ({\}) L [ ■ > 

On the other hand Theorem C.2 gives 

PJVU S > {A+}U{A_} 

Comparing (D.l) and (D.2) one can inductively deduce the form factors F^ N . The simplest way 
is to apply Theorem C.l which yields 

F i _ Np N ({\}) 

h L=0 

Let us introduce the matrix J which is given by pjy at L = 0: 



Jij = Sij ( ^2 ^ ik ) ~ 



k=i 

It is easy to see that 



L=0 i=l 



L 
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where the elements of J are the corresponding minors of J. The matrix J has the property that 
the sums of its rows and the sums of its columns give zero, therefore its minors coincide. Therefore 

FwAWn) = N 2 J u 

In graph theory the matrix J is the generalized Laplacian matrix of G n (the complete graph with 
n edges), where for each vertex from node i to j one associates the "weight" ipij. Then according 
to Kirchhoff's theorem the minor J\\ gives the enumeration of the spanning trees of G n : 

J\i = Gn(6i, ■ ■ ■ , Qn), 

where 

g N (e u ...,e N ) = \ U 

GGST \a€V G 

The sum runs over all spanning trees and the product runs over the vertices a of a given spanning 
tree, and 

<p(a) = tfij 

for a vertex a going from edge i to j. Putting everything together 

4w = n 2 ]r ( n 

GeST \aeVfc 

By the reverse application of Theorem 2.1 one can deduce that 

F 2N,c(^u ■■■An) = ^2 ^( A,J i ~ A <x 2 M A a 2 - A CT3 ) . . . ^(A CTJV _ 1 - A CTJV ) (D.3) 

Here the summation runs over all permutations a 6 Sjy. 

In the case of the operator O2 = W^f'^^B one makes use of the fact that it is just the 
interaction term in the Hamiltonian. Then the Hellmann-Feynman theorem gives [59] 

\ in L dc 

Here E(L, c) is the total energy of the iV-particle state for fixed momentum quantum numbers. It 
follows from scaling arguments that 

where h(x) is a dimensionless function. Therefore 

1 dE 1 fdE E 

= h 2— 

L dc c \ dL L 

The energy is given by 

^ ^ 9 , dE , d\i 

E = £ A? aid ^ = 2A,- 

i 

The derivative with respect to L can be calculated from the Bethe-equations: 



Ai + 0£ ~~ 
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Therefore 



\ 12 

\J/t\J/t^J/vJ> \ _ 

/ N pjv C L 



-■Jij + Sij ) Aj 



Making use of Theorem C.l one arrives at 



Fi N!S (e h ...,e N ) = -g N {6 u ...,e N )x\ E(a 4 - a,) 2 



Theorem 2.1 then yields 

-^VctAi,...^) = ^E^^i - A (T2 )v3(A (72 - X a3 ) . . .ip(X aN _ 1 - A CTJV )(A (71 -A CTJV ) 2 (D.4) 

The summation runs over all permutations a £ Sjy. The results (D.3) and (D.4) are in agreement 
with the formulas of [16] presented for N = 1,2, 3. 



E The x — > limit of the non-local correlation function 

In this appendix we demonstrate how to re-derive the expectation values {O^} from the pre- 
viously available results for non-local correlation functions [38]. We only consider O2 = ty^^ty , 
which is related to the operator Q 2 (x). The normal ordering of the field operators yields: 



px px 

Q 2 (x)-Q(x) = / dx 1 dx 2 ^(x 1 )^{x 2 )^{x 1 )^(x 2 ) 
Jo Jo 



Taking the limit x — > one gets: 



iJ/txT/txJnJ; _ lj m 



Q 2 {x) - Q{x) 



x- 



(E.l) 



Therefore, in order to obtain the expectation value of 2 = ^t^rt^r^r Q ne has to pick the 0(x 2 ) 
term from the expansion of Q 2 (x) — Q(x). We perform this for the expectation values with a finite 
number of particles. The relevant result reads (eq. 3.32 from [38]) 



N Pn({\)} 



[(Q(X)) 2 ) N + E ^({ A "}n, {\ + } N -n) P N-n({X + }) 
{A+}U{A_} 



(E.2) 



Here the summation runs over the partitions with |{A_}| = n with n = 2. . . N. The first term is 
given by the "irreducible part of the identity operator" (eq. 3.28 in [38]) 



((Q( X )) 2 )°N = E ^ Pn({K})p V N _ n ({\}) 



(E.3) 



{X x }U{Xy} 



Here /?^({A X }) and p y N _ n ({\ y }) are the Gaudin-determinants of a system with volume parameters 
x and y = L — x, respectively. The quantities / nj 7v({A_} n , {A+}Ar_ n ) are the "irreducible parts" 
of the operator Q 2 . They have a non-trivial dependence on x. In the case of Q{x} the "irreducible 
parts" are zero and the mean value can be expressed as (eq. 3.23 in [38]) 



<Q(*)> 



N 



Pn({\)} 



E n Pn({K})p V N - n ({X y }) 



(E.4) 



{A,}U{A y } 
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In the following we only consider the first non-trivial case N = 2. Equations (E.2), (E.3) and (E.4) 
yield 

Here 



(Q(x)Y - Q(x)) 2 = x 2pl (A l5 A 2 ) + 7 2 (Ai, A 2 ; 



p|(Ai,A 2 ) = x(x + 2ip 12 ) 

and 7 2 is given by 8.11-8.12 in [37]: 



2 



A 12 



( e -ix\ 12 _ ^ A12 + ic + ^^ Al2 _ ^ Ai 2 - ic 



Ai 2 - ic X12 + ic 

Expanding in x and keeping only the first two terms one has 



/ 2 = -4^ 12 + x 2 (^^-2)+... 
Putting everything together one observes that the 0{x) terms indeed cancel and one is left with 



{{Q{x)f-Q{ 



x 2 2lpi 2 \\ 



P2(Ai, A 2 ) 
Therefore 

12 p 2 (Ai,A 2 ) c 



Comparing with (3.17) we conclude that 

?2 _ n _ 2( ^12Ai 2 



F 2,s = °' ^4,s 



C 

The above formulas are in complete agreement with our result (D.4). 
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